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\Q . We present a parameter-free calculation of the dispersive and absorptive contribu- 

tions to the pion-deuteron scattering length based on chiral perturbation theory. 
We show that once all diagrams contributing to leading order to this process are 
included, their net effect provides a small correction to the real part of the pion- 
deuteron scattering length. At the same time the sizable imaginary part of the 
pion-deuteron scattering length is reproduced accurately. 

> 

X- 

1. The pion-nucleon (nN) scattering lengths are fundamental quantities of 
low-energy hadron physics since they test the QCD symmetries and the pat- 
tern of chiral symmetry breaking. As stressed by Weinberg long time ago, chi- 
ral symmetry suppresses the isoscalar irN scattering length a+ substantially 
compared to its isovector counterpart a_ . Thus, a precise determination of a + 
demands in general high accuracy experiments. 

Here pion-deuteron (jrd) scattering near threshold plays an exceptional role 
for Re(a 7r( i) = 2a + + (few-body corrections) . The first term ~ a + is simply 
generated from the impulse approximation (scattering off the proton and off 
the neutron) and is independent of the deuteron structure. Thus, if one is 
able to calculate the few-body corrections in a controlled way, nd scattering 
is a prime reaction to extract a + (most effectively in combination with pionic 
hydrogen measurements). In addition, already at threshold the nd scattering 
length is a complex-valued quantity. It is therefore also important to gain a 
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precise understanding of its imaginary part - this is one of the issues addressed 
in this letter. 

Recently the ird scattering length was measured to be [1] 

a^J p = (-26.1 ± 0.5 + £(6.3 ± 0.7)) x 10~ 3 m' 1 , (1) 

where m n denotes the mass of the charged pion. In the near future a new 
measurement with a projected total uncertainty of 0.5% for the real part and 
4% for the imaginary part of the scattering length will be performed at PSI [2]. 
Clearly, performing calculations up to this accuracy poses a challenge to theory 
that several groups recently took up [3,4,5,6,7,8]. In addition, an interesting 
isospin violating effect in pionic deuterium was found, see [9] . For a review on 
older work we refer to Ref. [10]. 

The imaginary part for the ird scattering length can be expressed by unitarity 
in terms of the ird total cross section through the optical theorem. One gets 

47rlm(a 7rd ) = lim q {a(ird -> NN) + a(ird -> 7 AA)} , (2) 

where q denotes the relative momentum of the initial ird pair. The ratio 
R = lim g _> ((y(ird -> NN)/a{ird -> jNN)) was measured to be 2.83 ± 0.04 
[11]. At low energies diagrams that lead to a sizable imaginary part of some 
amplitude are expected to also contribute significantly to its real part. Those 
contributions are called dispersive corrections. As a first estimate Bruckner 
speculated that the real and imaginary part of these contributions should be 
of the same order of magnitude [12]. This expectation was confirmed within 
Faddeev calculations in Refs. [13]. Given the high accuracy of the measurement 
and the size of the imaginary part of the scattering length, another critical 
look at this result is called for as already stressed in Refs. [14,15]. A consistent 
calculation is only possible within a well defined effective field theory — the 
first calculation of this kind is provided here. 

What is needed a priori for such an endeavor is a controlled power counting 
for NN — > NNir using chiral perturbation theory (ChPT) that is consistent 
with the one used for ird scattering. This was developed in recent years [16,17] 
- for a review we refer to Ref. [19]H This scheme led to the first calculation 
for pp — > dir + [20]. It was the central finding of this work that all loops 
that contribute to NN — > NNir at next-to-leading order cancel and the full 
transition amplitude up to next-to-leading order solely contains the diagrams 
shown in Fig. 1, however, with one modification compared to the standard 
treatment as used already in Ref. [21]: The irN — ► irN vertex in diagram ii) is 

1 For calculations that use chiral perturbation theory in the standard formulation 
for the reactions NN — ► NNir we refer to Refs. [18]. 
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i) ii) 

Fig. 1. Diagrams contributing to the ttNN — ► NN transition up to the order con- 
sidered: i) direct contribution and ii) rescattering. 

to be used with its on-shell value, 2m 7n instead of the previously used value 
of 3/2 m n . This increase was sufficient to bring the calculation in agreement 
with the data. 

In addition to the hadronic part of the dispersive and absorptive corrections 
to the nd scattering length, we estimate the corresponding contribution from 
the transition ird — > ^NN — > nd using the full structure of the one-photon 
exchange. Note that the inelastic channel ird — > ^yNN accounts for 1/4 of the 
imaginary part of a n d and therefore one can expect a sizable contribution also 
to its real part. 

The paper is organized as follows: in the next section we will present the power 
counting for the ird system including the dispersive part. In the third section 
we give our results, while a comparison to previous works is done in Sec. 4. 
We close with a brief summary. 

2. The basis of any effective field theory calculation is a proper power counting 
that allows one to organize the diagrams according to some a priori known 
hierarchy. It was argued by Weinberg [22] that in case of pion reactions on few- 
nucleon systems and especially on nuclei the Goldstone theorem ensures that 
one can expand the transition operators perturbatively. Then those have to be 
convoluted with proper wave functions. The important thing is to identify the 
relevant expansion parameter in the transition operators. For ird scattering in 
[22] the series is organized in powers of momenta and pion masses in units of 
the chiral symmetry breaking scale A x ~ 1 GeV. The typical one- and two- 
body diagrams are shown in Fig 2 (a) and (b), respectively. The small binding 
energy of the deuteron introduces a new small scale that can be accounted for 
systematically [3]. In Refs. [23,24] it was demonstrated how the scheme is to 
be modified in the presence of three-body (7rNN) cuts. Based on calculations 
with deuteron wave functions obtained solely from contact iViV interactions, 
in Refs. [25,26] it was argued that field theoretical consistency calls for a 
counter term at leading order. However, in Refs. [6,7,8] it was shown that this 
is no longer necessary as soon as the finite range of the one-pion exchange is 
included in the NN potential. 



a) b) c) 

Fig. 2. Typical Feynman diagrams for ird scattering; shown are the one-body term 
(a), the double-scattering contribution (b) and a four-nucleon-contact term (c). 
Solid black dots stand for the ttN interaction, whereas the hatched area shows the 
deuteron wave function. 

So far no attempt was made to also include consistently — within ChPT — the 
so called dispersive corrections that emerge from the hadronic nd — > iViV — > 
ird and photonic ird — > ^yNN — > red reaction chains. We define dispersive 
corrections as contributions from diagrams with an intermediate state that 
contains only nucleons, photons and at most real pions. Thus, all other poten- 
tially important corrections to the ird scattering length that come from, e.g., 
the A resonance will not be discussed in this letter. The diagrams contributing 
to the hadronic and photonic parts of the dispersive corrections in accordance 
to our definition are shown schematically in Figs. 3 and 4, respectively. 

Before we present the results of the calculation we first need to establish the 
power counting. The fact that the hadronic reaction chain ird — > NN — > ird 
is a process with large momentum transfer introduces a new scale into the 
problem that needs to be accounted for by a modified power counting. 

To establish the counting scheme we focus on two-body currents only - 
how to include one-body currents into the standard scheme is described in 
Ref. [220- Thus, in what follows we will compare our two body irNN — > 
NN — > irNN operators with the leading two body operator shown in Fig. 
2(b). Then it is sufficient to read off the vertex factors for the irNN — > irNN 
transitions to identify the order of any given diagram. We therefore estimate 
m 2 /(f^ q 2 ) for the diagram (b) of Fig. 2 where q here defines the momentum 
of intermediate pion. Utilizing Weinbergs counting scheme where all internal 
momenta are assumed to be of order m n we find diagram (b) to be (9(1) - 
here and in what follows we drop a factor 1//^ common to all diagrams to get 
the order estimate. Power counting gives that the 4N2ir contact term shown in 
Fig. 2 (c) contributes at 0(x 2 )> where x — m n/M N is the standard expansion 
parameter of ChPT with for the the nucleon mass. This last contribution 
comes with a yet unknown coefficient. As such, an estimate for its size provides 

2 In Ref. [16] the corresponding recipe is given that needs to be applied to NN — > 
NNir and therefore also to the dispersive corrections. It implies that the diagrams 
shown in Fig. 1(a) and (b) contribute at the same order for s-wave pion production. 
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Fig. 3. The classes of the hadronic contributions to the dispersive corrections. Di- 
agram (a) and (b) denote the direct and (c) and (d) the crossed terms. The filled 
ellipse denotes the NN interaction in the intermediate state. The diagrams con- 
tributing to the irNN — > NN transition to the given order are shown in Fig. 1. The 
diagrams with emission of pion on the second nucleon are not shown explicitly but 
taken into account in the calculation. 

the theoretical accuracy that a calculation for the ird scattering length can 
have at most. Therefore, assuming naturalness for the strength of the contact 
term, the theoretical limit of accuracy is of order x 2 which translates into a few 
percent. Reverting this statement, in order to reach a theoretical accuracy that 
is comparable to that expected for the experimental value of the nd scattering 
length, all contributions of lower order than x 2 should be evaluated. We will 
now show that the dispersive corrections contribute to 0(x 3 ^ 2 )- 

Transitions of the type ttNN —>■ NN —>■ ttNN — sketched in Fig. 3 (a) — get 
contributions from small values of the NN intermediate momentum q (q ~ m n 
or smaller) as well as from large values of q (q ~ p^ T = V m nM N ). The latter 
value refers to the on-shell momentum of the intermediate NN state (or, 
equivalently, to the threshold initial momentum of the reaction NN — > NNtt). 
The power counting as given for NN — > NNtt relates to the latter part of the 
contribution. It is based on the assignment p^/Mjy ~ 0(x 1 ^ 2 ) [16]. To derive 
the order of the dispersive corrections let us start with diagram dl of Table 
1. (See footnote 2 for how to include diagrams of the type d2). For this one 
we find in units of the amplitude for diagram (b) of Fig. 2 (estimated to be 

o(i)) 



jlqj m^~q 2 /M N + ie\{ATTf 




0( X 2 ) iorq 



O (x 2 I forg ~p thr 



(3) 



where the first term in the square brackets comes from the ttNN — > NN tran- 
sition operator, the second one from the NN propagator in the intermediate 
state and the last one from the integral measure. To arrive at the estimate for 
q ~ £>thr' wnere the l.h.s of Eq.(3) involves a singularity, we replaced the two- 
nucleon propagator by the corresponding 5-function term for this estimate 
should apply to both the real part as well as the imaginary part that emerges 



3 For a brief description on how to identify the order of a particular diagram we 
refer to Appendix E of Ref. [19]. 
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a) b) c) 

Fig. 4. Diagrams contributing to the dispersive corrections from photon-exchange 
interactions: one-body term (a) and double scattering (b). The filled ellipse denotes 
the NN interaction in the intermediate state. The solid squares denote the full 
irN — * jN transition amplitude as depicted in Fig. 5. 

from ttNN — > NN — ► ttNN. Furthermore we used 47rf n ~ Mat- The small 
momentum part of the integral is thus of order x 2 an d not relevant for this 
study. That is why dispersive corrections were not considered in the studies of 
Refs. [3,4]. However, the part of the integral where q is of the order of p^hr * s 
indeed of lower order than \ 2 an d thus should be considered. It is important to 
stress that for a consistent understanding of NN — > NNtt within ChPT it was 
also necessary to include the large scale p^r explicitly in the power counting 
[16,19]. For the imaginary part of the amplitude ttNN —>■ NN — > tcNN we 
have an experimental value — (3/4)Im(a 7r rf)/Re(a 7r rf) ~ 1/6, where the factor 
of 3/4 was introduced since this fraction of the width comes from ird — > NN. 
To check the power counting we need some estimate for the real part of the 
scattering length, which is known to be dominated by the double rescattering 
term (Fig. 2(b)) and was shown above to be 0(1). Therefore we expect from 
the above considerations a relative suppression of the imaginary part to the 
real part of the order of (m n /Mj^) 3 ^ 2 = 1/17. Thus, the hadronic contribution 
to the imaginary part of the nd scattering length is about a factor of 3 larger 
than predicted by the power counting — a deviation that is tolerable. 

Since the NN interaction is non-perturbative in diagrams, where a two- 
nucleon state contributes near on-shell, the full NN interaction should be 
taken into account — the corresponding diagram is depicted in Fig. 3(b). 
These are also included in our calculation using the techniques outlined in 
Ref . [24] . At the same chiral order there is an additional class of contributions 
- these are the crossed terms shown as diagrams (c) and (d) in Fig. 3. Al- 
ready in Ref. [27] an evaluation of those diagrams was called for, however, a 
consistent calculation of the terms that emerge from diagram if) of Fig. 1 has 
not been possible so far [10]. Since we work within a consistent field theory 
no such problem exists. The numerical importance of some crossed diagrams 
for ird scattering was already observed before and referred to as the so-called 
Jennings mechanism [28]: to understand data measured with tensor polarized 
deuterons for elastic nd scattering at backward angles, a particular crossed 
pion contribution needs to be included — see also the discussion in Ref. [29]. 
The order assignment for the crossed diagrams is obvious once one applies the 
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(a) (b) (c) (d) 

Fig. 5. Diagrams contributing to the ttN — > transition operator: Kroll-Ruder- 
mann term (a), pion in flight (b) and nucleon s- (c) and u- (d) channel. 

same procedure that leads to the estimate given in (3) — the only necessary 
change is to switch the sign of m n in the NN propagator. Note, in these dia- 
grams the two-nucleon intermediate state is always off-shell in contrast to the 
NN state in the direct contributions that are expected to receive the dom- 
inant contributions from (near) on-shell nucleons. It is therefore surprising 
that direct and crossed terms appear at the same order. However, one should 
recall that the chiral expansion is an expansion around the chiral limit. When 
approaching the chiral limit the intermediate two-nucleon states of both di- 
rect and crossed diagrams approach the same kinematical point. Therefore it 
is natural that they also contribute to the same chiral order for physical pion 
masses. 

A priori there is no rule how to include the electromagnetic contribution to 
both the real and the imaginary part of the ird scattering length (see Fig. 4) 
into the power counting — the fine structure constant a is clearly an indepen- 
dent parameter. Based on the observation that the electromagnetic and the 
hadronic contribution to the imaginary part are of the same order of magni- 
tude, we assign the same chiral order to both — as it is often done in chiral 
perturbation theory studies. 

3. In this section we present the results of our investigations. We first focus 
on the hadronic contribution to the dispersive corrections. The results are 
given in Table 1. Note, for the contributions that involve the NN interaction 
in the intermediate state we do not give the individual results explicitly. All 
calculations are done with the CD-Bonn NN potential [30]. 

First of all we observe that with a value of 4.25 x 1CT 3 m" 1 the imaginary part 
of the scattering length turns out to be very close to the experimental number 
of (3/4)Im(a^ P ) = (4.7±0.5) x 1(T 3 m' 1 (c.f. Eqs. (1) and (2) and discussion 
below the latter). This is not too surprising given the good description for 
the near threshold cross section of NN — > dn as reported in Ref. [20]. For 
the real part on the other hand we observe a striking pattern: although the 
individual contributions can be quite large, the total sum turns out to be rather 
small. About 1/3 of the contribution from the direct diagrams of Table 1 gets 
canceled by the corresponding ones with NN interactions in the intermediate 
state. This is the same pattern as for the imaginary part — the fact that the 
latter reduction in the magnitude is natural for processes of the type NN — > 
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NNx was discussed in Ref . [31] . However, about 60 % of the contribution of the 
direct diagrams of Table 1 is canceled by the corresponding crossed diagrams. 
As we will see, part of this cancellation is quite natural. When comparing the 
direct and the crossed diagrams we observe that the two-nucleon propagator in 
the intermediate state of the former reads l/(m n — q 2 /M^) (c.f. Eq. (3)), where 
q denotes the relevant loop momentum. The corresponding expression for the 
latter reads l/(— m n — q 2 jM^\ Thus, for small values of q, where the two- 
nucleon propagator becomes static, one obtains l/m n and respectively, 
and some contributions of direct and crossed diagrams, specifically d2 and 
c2, will largely cancel. Note, this cancellation does not mean that the full 
contribution of each pair of diagrams cancels. The numbers given in the table 
also contain the contributions from large values of q, where such a cancellation 
does not necessarily occur. Note also that it is the structure of the two nucleon 
propagator that is responsible for the smallness of diagram dl as compared to 
its crossed partner cl. In contrast to diagram cl, the two-nucleon propagator 
in dl changes its sign when passing through the point q 2 = m^M^ whereas 
all other terms such as vertex functions, pion propagators, etc., have the same 
sign throughout the region of integration. Thus, a strong cancellation takes 
place for the latter and as a result the real part of diagram dl is much smaller 
than that of cl. The situation for diagram d3 is different to dl, since also the 
S-wave deuteron wave function changes sign at q ~ y/m^M^. This leads to a 
constructive interference of contributions from small and large momenta. 

In Fig. 4 we show the diagrams that contribute to the electromagnetic piece of 
the dispersive contributions. To evaluate the real part of the one-body contri- 
bution (diagram (a)) we use the same prescription as used in Ref. [23], namely 
we subtract the term corresponding to the one-body operator at zero momen- 
tum. This removes the leading divergence that in a full calculation needs to 
be absorbed into the electromagnetic corrections of the irN scattering lengths 
(note, this quantity was recently estimated in Ref. [32]j3. To ensure gauge in- 
variance the ttN — > 7iV amplitude needs to contain all diagrams shown in Fig. 
5. However, since for the calculation of the photonic part of the dispersive cor- 
rection the tcN — ► 7iV amplitude contributes at threshold, diagrams (c) and 
(d) are suppressed by one power of x as compared to diagrams (a) and (b). 
Thus, we did not take them into account in this leading order calculation. Be- 
cause of the same reasoning we neglect also the contributions to the irN —>■ 'jN 
amplitude from the A resonance [34] . We calculated the full red —>■ '-/NN — ► ird 
amplitude in Coulomb gauge. The corresponding propagator is given, e.g., in 
Ref. [3 5]. The final result is (—0.1 + i 1.4) x 10~ 3 m^ 1 . Again, the imaginary 
parlll is sufficiently close to the corresponding experimental value of 



A more precise calculation within QCD+QED requires a much more sophisti- 
cated framework, see e.g. [33] — this is beyond the scope of this letter. 
5 Note that the photonic part of the absorption correction is dominated by the 
imaginary part of the 7r~p scattering length due to the process ir~p — > — > ir~p. 
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Table 1 

Hadronic contribution to the real and imaginary part of a^d m units of va~ x x 10 -3 . 
We only show the typical topologies — all permutations are included as indicated 
by the ellipses. Line code as in Fig. 2. 



direct contributions 



dl 



d,2 



d3 



= 



N + 4 





r + 



r + 



0.49 + i6.68 



-0.8 



-6.02 



i 0.1 



i 1.65 



dA-dG 



sum of this group 

corresponding terms with intermediate NN interaction 



-6.33 + ib.13 
1.96 - i0.88 



sum of all direct terms 



-4.37 + i4.25 



crossed contributions 



cl 



c2 



c3 






1.83 



1.37 



-0.09 



c4-c6 



sum of this group = 3.11 

corresponding terms with intermediate NN interaction = —0.37 



sum of all crossed terms 



2.74 



direct + crossed 



total sum 



-1.63 + i4.25 



9 



(l/4)Im( a nd J = (1.6±0.2)xl0 3 m^ 1 , whereas the real part is negligible. 
This conclusion is consistent with that of Ref. [14], where the real part was 
found to be —0.2 x 10~ 3 m~ l within a less complete calculation^!. 

4. We now compare our results for the hadronic contributions to other works. 
The dispersive corrections to ird scattering were investigated using Faddeev 
calculations in Refs. [13]. Since these works considered only those intermedi- 
ate states that contain at most one pion at a time, all direct diagrams were 
included. The authors found —5.6(1.4) x 10 -3 m" 1 as contribution to the real 
part — this number was also used in the reanalysis of ird scattering in Ref. 
[14]. This is to be compared with —4.37 x lO^m^ 1 from our work, which 
agrees to the above result within the uncertainty. Note, the Faddeev equa- 
tions produce amplitudes that are non-perturbative in the irN and the NN 
interaction simultaneously. Thus, in addition to the direct terms as shown in 
Table 1 there are many more diagrams included like contributions where the 
iViV-pair interacts while there is a pion in flight. However, all those are of 
higher order in the chiral expansion. The closeness of our result for the direct 
terms and the corresponding result of Ref. [13] is thus an indication for the 
convergence of the chiral expansion. Recall that our final result to the real 
part of a-nd is smaller because of a cancellation of the mentioned contribution 
with the crossed diagrams that were not included in Refs. [13]. 

In Ref. [5] the diagrams dl of Table 1 were evaluated explicitly, besides many 
others that are difficult to match to our amplitudes. In this work they con- 
tribute with (0.24 + 26. 0) x lO^m" 1 to the ird scattering length. This number 
is to be compared to our (0.49 + z6.68) x lO^m' 1 . In Ref. [5] a value of 
3/2m n was used at the ttN — > tcN vertex in contrast to the proper 2m n as 
derived recently [20]. On the other hand, Ref. [5] includes a (small) isoscalar 
piece into this vertex — chiral power counting assigns a subleading order to 
this piece. However, the pattern of the result is the same: a small real part is 
accompanied by a sizable imaginary part. 

Our results for diagrams d2 and c2 agree to those of Ref. [37], once the same 
ttNN coupling constant is used. 

Some of the diagrams in Table 1 where included in the phenomenological stud- 
ies of Refs. [5,15]. In particular, the second diagrams of d3 and c3 contribute 
to the so-called SP-interference term (a double scattering diagram, where the 
first ttN — > ttN transition is in an s-wave, whereas the second is p-wave). 
In those studies the p-wave amplitude was taken from fits to irN data and 
parameterized as a strength parameter times the square of the irN relative 



6 In Ref. [36] the corresponding integral was evaluated to be one order of magnitude 
larger; however, the author assumes the photon energy go in the exchange to vanish 
instead of using the correct value of qo ^ m w . 
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momentum. Indeed, what appears at the lower nucleon can be regarded as 
the 7riV scattering in £>-wave but in the boosted frame. However, this treat- 
ment misses an important momentum dependence, since in the boosted frame 
the nucleon propagator in the intermediate state contains a term q 2 /M^ in 
the denominator (see Appendix). As a consequence the p-wave subamplitude 
in the phenomenological studies grows quadratically with momentum, even 
for momenta q 2 ~ m^M^ and the full matrix element scales with the wave 
function at the origin [5], which is theoretically not controlled. On the other 
hand, in our case this subamplitude goes to a constant leading to a controlled 
behaviour of the matrix element [7]. This is why those parts of Refs. [5,15] 
can not be matched to our results. Based on the arguments given, we call for 
microscopic calculations, now possible within ChPT, instead of applying the 
phenomenological techniques. 

5. To summarize, in this work we calculated for the first time the absorptive 
and dispersive corrections to the nd scattering length using ChPT. Especially 
we found for the absorptive part 

Im(a^) = ((4.25 ± 1.2) + (1.4 ± 0.4)) x 1(T 3 m~ l (4) 
to be compared with the experimental value of 

Im (a^ P ) = ((4.7 ± 0.5) + (1.6 ± 0.2)) x 10~ 3 m' 1 . (5) 

In both Eq. (4) as well as Eq. (5) we give the hadronic and electromagnetic 
contribution separately. We thus find good agreement between theory and ex- 
periment for each of the contributions. The theoretical uncertainty is estimated 
to be of order 2m 1T /M N in both cases[3. 

For the corresponding dispersive part we get 

47 = -1-7 x lO- 3 ™; 1 =► off/Re («^ P ) ~ 6.5% . (6) 

The number given for a^ p now contains both the hadronic as well as the 
electromagnetic contribution and for Re(a 7rd ) we used Eq. (1). This result 
is quite small given that the imaginary part of the ird scattering length is 
about 1/4 of the real part. It is difficult to provide a proper estimate for 
the uncertainty of a^ p \ since it emerged from a cancellation of individually 
sizable terms. The most naive method would be to use the uncertainty of 
order 2m 7r /Mjv one has for, e.g., the sum of all direct diagrams to derive a 

7 The factor of 2 appears, because the ttNN — > NN and irNN — > ^yNN transition 
operators — both evaluated with an uncertainty of order m n /M]y — appear twice 
in each amplitude. 
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Aa^ sp of around 1.4 x 1CT 3 m^ 1 , which corresponds to about 6% of Re (a^P) • 
However, given that the operators that contribute to both direct and crossed 
diagrams are almost the same (see Appendix) and that part of the mentioned 
cancellations is a direct consequence of kinematics, this number for Aa^f 
is probably too large. A reliable error estimate for a^ p requires an explicit 
evaluation of the next order contributions. 

We showed that the smallness of the dispersive contribution to the real part of 
the nd scattering length is a consequence of efficient cancellations amongst var- 
ious, individually sizable terms. To gain a better understanding of the real part 
of the pion-deuteron scattering length, a complete calculation of all isospin- 
breaking corrections at N(N)LO in the EFT with virtual photons is called for 
(as also stressed in [9]). 
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Appendix 



In this appendix we present the explicit expressions for the amplitudes de- 
picted in Table 1. Note, in accordance with the definition used for dispersive 
corrections as well as the power counting, we only keep those amplitudes that 
contain two-nucleon cuts in time-ordered perturbation theory (TOPT). Es- 
pecially, we dropped the so-called streched boxes. 

Using the same labels as in the table, one finds for the corresponding correc- 
tions to the nd scattering length 



dis P = g\ m l f d 3 q q 2 (a+[3) 2 _ , , 

nd 6nf«(l+m n /2M N )(l+q /2M N ) J (2ir)* q - q>/M N + ie 1 " ' 



where 
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a =(/,(,)- ^/,( 9 )) and = -f («(,) + . (A-2) 

Here ii and .Z2 are the integrals that correspond to the overlap of the deuteron 
wave function (u(q) and w(q) for the S- and D- waves, respectively) with the 
one pion exchange operator 




'*>W|fr (1 - tf ff W)) -< + *> 

where Pi and P2 correspond to the TOPT components of the pion propagator 
Pi = g - - (p 2 + q 2 )/2M N and P 2 = -uJ^+q - (p 2 + q 2 )/2M N with 
= ^q 2 + m\. 

The diagrams of Table 1 can be easily matched to the individual terms of 
Eq. (A-l): type 1 contains a 2 , type 2 contains f3 2 , whereas the interference 
terms of type 3 contain 2a(3. For the direct terms, labled as d in the table, 
one needs to take g = m w and for the crossed terms, labled as c in the table, 
go = —m n . All necessary information on how the AW interaction is included 
in the intermediate state can be found in the Appendix of Ref. [24]. 

All integrals are evaluated up to a sharp momentum cut-off of 1 GeV. All 
higher momentum components are to be absorbed in a counter term that is to 
be included at order \ 2 ( c -f- discussion in section 2). By enlarging the cut-off 
by a factor of 20 we checked that the integrals change by less than 10 % - 
fully in line with the power counting. 
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